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Abstract 



We study the production and the equihbration of a non-Abehan qq plasma 
in an external chromoelectric field, by solving the Boltzmann equation with 
the non-Abelian features explicitly incorporated. We consider the gauge group 
SU{2) and show that the colour degree of freedom has a major and dominant 
role in the dynamics of the system. It is seen that the assumption of the so 
called Abelian dominance is not justified. Finally, it is also shown that many 
of the features of microscopic studies of the system appear naturally in our 
studies as well. 
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I. INTRODUCTION 



There is a high expectation that with the advent of new machines (RHIC at BNL, LHC 
at CERN) it will be possible to study the quark gluon phase of strongly interacting matter. 
Lattice QCD predicts that at density and temperature of the order of 2GeV / frn? and 200 
MeV respectively, hadronic matter undergoes a transition to the quark-gluon phase the 
Au-Au collisions at 200 GeV/nucleon (in the centre of mass frame) in RHIC are expected to 
provide such energy densities which would then give rise to the quark- gluon phase [Q]. It is 
also generally believed that the transition is of the first order, although more recent lattice 
studies IP] and those of intermittency suggest that this belief needs more substantiation. 

The most important question experimentally is of course the signature of quark gluon 
plasma (QGP). Many proposals have been made, prominent among them being J/\E' suppres- 
sion strangeness enhancement [§] and dilepton production [0. Clearly, it is important 
to know, even if only roughly, the dynamics of the QGP evolution in space - time. Indeed, 
apart from being intrinsically interesting, the space - time dynamics will have a direct impact 
on the signatures proposed. 

The space time evolution of QGP that is produced in ultra-relativistic heavy ion 
(URHIC) collisions can be described P broadly in four stages. They are, i) pre-equilibrium, 
ii)equilibrium, iii) cooling, and iv) hadronisation. While most of the earlier work refers to 
stage (ii) where one studies the hydrodynamic evolution 0, subsequently there has been a 
good deal of effort in understanding stage (i) as well. Since it is of relevance to us here, we 
briefly reveiw them below. 

In the context of Flux-tube model (described in the next section), the transport equation 
is solved to study the evolution of the system. Assuming the Bjorken scenario [|10l where 



a central rapidity plateau will be seen in relativistic heavy ion collisions |]rT|, this model 
deals with a baryon free plasma in the central region; there would initially be a huge deposit 
of energy which then creates the partons. In the approach of Baym these partons 
are assumed to be created soon after the two nuclei have collided with each other, and are 



further assumed to equlibrate almost instantaneously thereafter (with a very small relaxation 
time). This quark-gluon gas then proceeds to expand hydrodynamically. Subsequent work 
of Kajantie and Matsui [0], incorporated a dynamic particle production by introducing a 
source term via the Schwinger mechanism |T^ in the Boltzmann equation, and they also 
studied the equilbration in the relaxation time approximation. On the other hand, Bialas 
et al 1 14 1 generalized the Baym analysis by introducing the effect of the background field on 



the otherwise hydrodjTiamic flow. A self consistent study of the system was carried out by 
Banerjee et al who combined the effect of both the background field and the collisions on 



the quark-antiquark system. Note that in the analyses of Kajantie et al |12] and Banerjee 
et al [0], the source term acquires a time dependence by virtue of the time dependence 
that the electric field suffers because of the particle production. Finally, by employing the 
same analysis as of Banerjee et al, Asakawa and Matsui [|16[ have studied the rate of dilepton 
production by qq annnihilation, as a function of (proper) time. All the above analyses ignore 
the gluon component. 

In contrast to the above studies which we shall broadly term macroscopic, there is an 
interesting and a more complete as well as a complementery approach by Geiger and Ka- 
pusta [T^. They have studied the QGP evolution in what they call the parton cascade 
model. In this approach, the production and the evolution of qq pairs and gluons are stud- 
ied numerically, by considering the direct collision processes 2 — > 2 , the inelastic fusion 
processes 2^1, and the decay processes 1-^2. The above processes are studied within the 
frame work of the perturbative QCD, by taking the leading order and the next to leading 
order contributions. While this approach of Geiger and Kapusta has the merit of being 
microscopic, it comes with the complication of extensive numerical computation - and the 
physical picture may not always be clear. More significantly for our purposes, the analysis 
confirms the Bjorken hypothesis of the boost invariance of the parton densities in the central 
rapidity region. Further, the relaxation time which is other-wise free in the flux-tube models 



also gets determined. Indeed, Geiger et al [|I^] report a value of ~ 0.2/m in their paper. 
One may thus look upon the flux tube models to be effective versions of the more detailed 



microscopic pcm, if the Tc determined by the latter is employed in the transport equation. 
We shall make a more detailed comparison in the subsequent sections. 

In any case it is our purpose to emphasize and study a crucial feature of the quark- 
gluon phase, viz, the inherently non- Abelian nature of the plasma. Note that all the 
earlier macroscopic approaches ||^,|T^,|T^-|TC[ are in the so called Abelian approximation. Put 
simply, the non-Abelian features are completely ignored in the space time evolution, and 
the Boltzmann equation describes the dynamics which is essentially that of a Maxwell fluid, 
be it in the source term or in the background field term. The approach of Geiger and 



Kapusta [17| does incorporate the non-Abelian features in their basic Feynman diagrams, 
but is neverthless incomplete in that the distribution function in its final form does not carry 
any color degee of freedom. 

We attempt to fill this gap, within the framework of transport equations, by explicitly 
incorporating the non-Abelian features mentioned above. To be sure, the relevant formalism 



has already been devloped by Heinz [18| . Here one needs to carefully implement the non- 
Abelian constraints within the Boltzmann equation such that the Bjorken picture of boost 
invariance along the collision direction is also preserved. 

Before we proceed to present the formalism it is pertinent to make the following obser- 
vations: 1) Strictly speaking, the source term that has been employed in the above studies 
is inadmissiible. For, even within the so called Abelian approximation, the chromoelectric 
field suffers a space time evolution in any self-consistent calculation. To wit, the production 
of particles is then strictly governed by a time dependent field, in which case a perturbative 
mechanism will take over from the non-perturbative Schwinger mechanism. Note that the 
latter holds only for a uniform constant field. All the previous studies employing the flux 
tube model employ the non-perturbative expression. 2) As already mentioned, the source 
term is derived for QED processes and needs to be redone for QCD processes at hand. It is 
apparent that the QCD effect will show up most manifestly in the gluon production; there 
is no corresponding counterpart in QED. In other words it is necessary to obtain a source 
term for gluons. 3) It is well known from earlier studies of Yang-Mills equations that there 



are inequivalent gauge field configurations which yield the same field tensor ||I9[, and hence 
the same energy momentum tensor. This feature which is again absent in the Maxwell 
case has to be taken into account in our study. Recall that the lattice studies at any given 
energy density implicitly sum over the contributions coming from all such inequivalent con- 
figurations 4) Finally, the inherent non-linearity and the local gauge invariance which 
characterize the Yang-Mills equations lead to an added complexity; the equations do not 
possess a unique solution, even after gauge fixing, and even if one considers only static con- 
figurations. Consequently, one needs to scan, with some suitable weightage, all the solutions 
for the given initial conditions, and sum over all such configurations. It is not clear if any 
procedure is known that allows an implementation of this requirement. 

In this first study we restrict ourselves to a study of the simplest of the situations leaving 
the points mentioned above for a later presentation. Apart from acting as a warm up, it will 
also allow us to devise a simple bench mark that would enable us to compare the Abelian and 
the non- Abelian features. With that in mind we consider only Maxwell like fields (see below), 
ignore the production due to a time dependent field and, finally restrict ourselves to the gauge 
group SU{2). The last choice leads to an additional simplification as quarks, antiquarks and 
gluons all belong to the same represantation. As a final simplification, we shall also drop 
the gluonic terms altogether. Please note that even with all these approximations (which 
we shall remedy in later publication) our approach is none-the worse compared to the fully 
Abelian models which have been employed so far. 

The plan of the paper is as follows. In the next section, we shall briefly describe the 
model where we set up the transport equation with the non-Abelian features appropriate 
to the group SU{2) explictly incorporated. We also obtain the auxiliary equation that 
follows from the requirement of the conservation of the energy momentum tensor. Section 
3 describes briefly the numerical procedure employed. In Section 4 we present the results 
and discussions. We conclude and summarize the main results in Section 5. 



II. THE FLUX TUBE MODEL 



One well known signature for the QGP is the dilepton production which has the 
attractive feature that the leptons do not suffer any final state (strong) interaction. As 
such this process can therefore be used to study the space-time evolution of the QGP. 
In particular, one would be able to study the pre-equilibrium regime, provided one can 
determine the distribution functions for the quarks and the anti-quarks at such early times. 

The flux tube model, as we have mentioned, is a useful choice in this context. We shall 
employ this model in this paper. A brief description of the Flux Tube model is not out of 
place here. 

In this model the two nuclei that undergo a central collision at ultra high energies are 
lorentz contracted as thin plates [pl| , p2| . These two lorentz contracted nuclei pass through 



each other and, in the process, acquire a non-zero colour charge {{Q) = 0, (Q^) ^ 0), by 
exchanging soft gluons. So one may figuratively call such nuclei after collision as color 
capacitor plates between which a strong chromoelectric field is created. This field creates qq 
and gluon pairs via the Schwinger mechanism which enforces the instability of the vacuum 
in the presence of an external electric field. The partons so produced will collide with each 
other and also get accelerated by the parent background field. The mutual collisions drive 
the system towards equilibrium with suitable modulation from the background acceleration. 

Let /(/) be the one particle(anti-particle) distribution in the phase space. In the earlier 
Abelian approximation, the evolution of the particles is studied via the Boltzmann-Vlasow 
equation, 

},d^ ± gF.^p'^d^] fix,p)ifix,p)) = Cix,p) + S{x,p), (1) 

where g is the strength of the charge which couples to the field. Note that the distribution 
functions /, / are defined in a six dimensional phase space of coordinate and momenta. In 
the above equations, C is the collision term driving the system into equilibrium, and S is 
the term that acts as the source for pair production |T^. We now generalize Eqn. (1) to 



incorporate the non-Abelian features explicitly. An easy way to to accomplish this is to 
extend the phase-space by taking the color degrees of freedom into account. The extended 
phase space is now a direct sum IZ^ © Q, where Q is the (compact) space corresponding to 
the given gauge group. Thus the single particle phase space has a dimension d = 6 + (A^^ — 1) 
if we consider SU{N). The evolution of the color charge now follows Wong's equation 

^ = r^'u.Q'A^'^ (2) 

which supplements the Lorenz force equation 

P ^ Qa;pa^.^^ (3) 

where is the four-velocity and f"''"^ are the structure constants of the gauge group. 
For SU(2) charges that we are interested in, the non-Abelian extension of Eqn.(l) reads 

p,d^ + Q'F^^yd^ + e^'V^'Ayd'ol fix, p, Q) = C{x, p, Q) + S{x, p, Q) (4) 

Note that we do not have to write a separate equation for anti-quarks since they belong 
to the same color representation as quarks. Indeed, the antipodal points on the sphere 
(corresponding to the color part of the phase-space with a fixed charge) represent the particle 
and anti-particle. Eqn. (3) may be recast in the more convenient form 

^=u,Qx (5) 

where the arrows now denote the direction in the color space. In the same notation, the 
transport equation obtains the form 

d 



f{x,p,Q) = C{x,p,Q) + S{x,p,Q) (6) 



dQ\ 

For the model at hand, let us take the 'plates' to be moving along the z-direction. As 
mentioned earlier, we restrict F^^ to be Maxwell like, i.e., we restrict to only that field 
configuration for which there exists a gauge choice such that the gauge potentials commute 
with each other every where. We next require a boost invariant description of the distribution 



functions as well as the other physical quantities which may be determined thereof, even 
as the system is evolving. Keeping this in mind, we make the gauge choice where only 
the components A'*" = [A^^,A^^) are non- vanishing and the other components are zero. 
With this choice, the resulting electric field points in the '3' direction in the colour space. 
We now proceed to impose the Lorcntz gauge condition, which has to be done such that 
the chromoelectric field depends on the boost- invariant (along the axis of collision) quantity 
T = {t^ — z^y^'^. Observing that the system is effectively (1+1) dimensional, in the t-z plane, 
we write 

A^" = e^^a^G(r) (7) 

where the indices /i., v take values 0, 3 and G is a Lorentz scalar function depending only 
on T. The above choice automatically implements the Lorentz gauge condition d/^A^^ — 0. 

— * 

Clearly, the chromoelectric field E is dependent only on r, and is given by 



E(r) 



dr^ ^'/^\r 



G{r) (8) 



The Wong equation guarantees the conservation of the magnitude of the vector charge 

— * 

Q, which may now be held fixed. Being the analogue of the Larmor equation for a charged 
particle in an external magnetic field, it also conserves the component of the charge that is 
parallel to the external chromoelectric field. It is therefore convenient to resolve the SU(2) 
charge in the polar coordinates. Writing 

Qi = QsinOcoscf); Qi = Qsin9sin(j); = QcosO , (9) 

it is straight forward to verify that 

The last equation will lead to considerable simpification in solving the transport equation. 
Let us now consider the collision term. As mentioned in the introduction, the colour 
plate model can be looked upon as an effective version of the more detailed parton cascade 



model which has confirmed the Bjorken hypothesis. We may therefore, employ the simple 
relaxation time hypothesis, with the phenomenological parameter Tc, the relaxation time, to 
be obtained from microscopic computations. It then follows that 

^ _ ~P^'^^lif ~ /eg) {ii) 

where feq is the eqilibrium distribution function, with local (space time dependent)values 
of the thermodynamic quantities. While it has been customary to take feq to be one of an 



ideal gas, with some support from the pcm analysis of Geiger and Kapusta ||T^, recent lattice 
studies suggest that the quark-gluon phase is possibly a non-perturbative phase. Since we 
are not dealing with a true system here, we shall conveniently take feq to be that of an ideal 
Fermi gas with a local temperature, evolving as a function of the proper time r. We thus 
have 

2 

^""^ ^ exp(;j/^V^(r)) + l' ^^^^ 
where is the flow velocity, 

= {coshrj, 0, 0, sinhrj). (13) 

which is written in terms of the (space-time) rapidity tanhrj = z/t. 

Now we demand the boost invariance following Bjorken's picture according to which the 
longitudinal boosts are the symmetry operations on the single particle distribution. The 
boost invariant parameters on which / can depend are, apart from the charge coordinates, 

r=it'- z^f/\ ^ = iv-y),Pt = {pl- V]f" (14) 

where y = tanh~^(pi/]9o) is the momentum rapidity. It is convenient to write separate 
equations for quarks and anti-quarks, a la the Abelian case. If we therefore identify the quark 
states with the points on the upper hemisphere of the colour sphere, and the antiquarks with 
the lower hemisphere, by a trivial relabeling, we may write two equations. 



d /tanh,^ gcoi^OEir) 
dr \ T ptcosh^ 



- + ,-G(r)tanhe^ 



+ 



eq 



Ptcosh^ 



+ 



d /tanh^ gcos6E{T)\ d ^ ^ G( )t h^*^ 
dr \ r ptcosh^ J 8$, dr d(j) 



Tc Tc Pt cosh ^ 



(15) 
(16) 

(17) 
(18) 



with the first of them for the quarks and the next for the anti-quarks. The angle variable 
9 varies from to 7r/2 in both the equations. Finally, S is the (non perturbative) Schwinger's 
expression for pair production and is given by 



V/ c Q\ gEcosO 
^{r,tPt,9)^ ^^^^'^ 



1 — exp 



TT 



(19) 



gfE'cos^^ 

where we have inserted the Gaussian dependence on ^ by hand. Note that the non- 
Abelian nature features dominantly via the cos{9) term. This occurence plays an important 
role in the evolution of the system. 

After having set up the relevant equations (not all yet since energy-momentum conserva- 
tion is to be imposed), we observe that the above differential equation possesses the formal 
solution 



j-T q-' — q- 

f{r,^,Pt9,(p) = o?T'exp( ) 

Jo Tc 

where Uri) is given by 

i' = sinh"^ 



^{r',C',Pt,9) ^ feq{r',C',Pt) 
Pt cosh ^' Tc 



(20) 



r , ^ qcos9 r , ,, ,,, 
- sinh i + / dT"E(T") 

t' PtT' Jt' 



(21) 



It is thus clear that / does not depend on 0, so the ^ term contributes nothing to the 
transport equation. However the 9 dependence is still involved through out the formalism 
which shows the non-Abelian effects. The corresponding distribution / for antiquark can be 
written by just changing g to —g in / as can be checked out easily. 

The transport equation has to be supplemented with the conservation constraint 



dfiTlf^at + (^plTym — 0) (22) 

since it is the field energy that is being pumped in order to produce the qq pairs. More 
exphcitly, 

d.TZt = -r.Fr ^ -d.Tpi, (23) 

with 

TZt = jv'fU + DdVdnQ (24) 

where the measures dVtn = smOdedcj) , dV = = and = J p^'Qaif - 

f)dTdVLQ. Here we have taken the value of the degeneracy factor 7 = 2 corresponding to 
two flavours. Since energy and momentum are conserved in each collision, the moment of 
the sum of the collision term vanishes: 

p''{C)dVdVLQ = (25) 

Now taking the first moments of the Boltzmann equation and integrating over the color 
degrees of freedom for / and / and making use of the conservation of energy and momentum, 
we obtain 



a^Tf + gE{r) J dTdQQp'^ ^^-^ + 2 J dTdQQS = (26) 

where 

Tj^" = diag(EV2, ^V2, ^V2, -£^2) (27) 

is the energy momentum tensor for the field. We may solve for the electric field by 
employing the same procedure as in |jl5|, and by employing the symmetry f{T,C,,PT,0) = 



f{T, —C,,Pt, 6)- We thus obtain the equation governing the decay of the source field to be 

^ - rdptpl r^esinhe r"de,m6[f - /] + ^-a\E{r)f' = (28) 
ar Ztt Jo Jo Jo t 

where a = aC(5/2) exp(0.25/«), a = c{g/2f'^j^ and c = Finally, C(5/2) = 1.342 

is the Reimann zeta function. 



Equations (37) and (26) are as yet underdetermined since the local temperature T(r) is 
free. In order to fix the form of T{t), we appeal to the relaxation time approach that we 
are employing, and assume that the the particle energy density differs negligibly from the 
equilibrium energy density. We may then relate, by an ansatz, the proper energy density 
which is defined by 

e(r) = JdTdnQ{p^u,y{f + f) (29) 
to the temperature by its equilibrium value, whence. 

It may be mentioned that the weaker condition of energy- momentum conservation that 
we have employed here , in fact, satisfies the Yang-Mills equations D^F^^''^ = j" as well. 

Finally, we pause briefly to discuss the effects of hard thermal loops that have been 
emphasized recently [Q. They have been derived from the classical transport equation as 
well by Kelly et al P3| , |2B[] . The latter derivation, which is important for us here, is based 
on an analysis of the Vlasov equation which properly describes the expansion of an already 
equilibrated gas, in a back-ground field but with with no source or sink. It is, therefore, 
necessary to determine whether the more complete transport equation such as the one that 
we have here will have any effect on the results of Ref. ||25|j2^. Conversely, it is also necessary 



to determine how the hard thermal loops will affect the purely classical results that will be 
obtained here. Let us recall that the derivation of Kelly et al |2^j2^ consists of a systematic 
expansion of the distribution function as well as the background field term in powers of the 
coupling constant. The zeroth order term is merely the Fermi-Dirac term for the quarks. 
Since we are solving here the approach to the equlibrium, there is necessarily a dependence 
of the distribution functions and other physical quantities on g. It remains to disentangle 
the contribution coming from the hard thermal loops. Indeed, observe that the two extra 
terms that we have at hand here are of higher order in g. First of all, the source term is 



non-analytic and does not even admit an expansion in powers of g The other collision 



term is easily seen to be of order g"^ or higher. We thus conclude that the conclusions of 
Ref. |]25|j2^ remain unaffected, and we may take over those results in toto, supplementing 
our classical results. 



III. COMPUTATIONAL PROCEDURE 

We have adopted here a double-self consistent method to determine (/, e, ...) , following 
the work of Banerjee et al [|l5l in the Abelian case. The procedure follows the scheme 
{TiT)triah E{T)triai} ^ {/, / E{t)} {/, /} ^ T(r) ^ ... by repeated use of equations 
(20), (28), (20), (30). The iteration terminates as soon as a convergence is established in 
the solutions for E{t),T{t). All the desired quantities are thereby consistently determined. 



IV. RESULTS AND DISCUSSIONS 

Before we proceed to present and discuss the results, a few comments about the choice of 
the value of the parameters in the computation. We put g = 4: throughout our calculations. 
Since lattice computation results predict a phase transition from the baryonic phase to the 
quark-gluon phase at densities ~ 5 — 10 GeV, it has also been customary to take an initial 
energy density in the same range in the Flux tube model. However, greater care needs to 
be taken before the initial energy densities in URHIC are chosen. Indeed, a fairly reliable 
estimate of the time required for the partons to be produced in central collisions after the 
two nuclei have suffered the maximum overlap is ~ 0.05 — 0.1 /m; a simple dimensional 
analysis leads to a value of the initial energy density to be {1/2)Eq ^ 500 — 1000 GeV/ fm^, 
the precise value depending on the the magnitude of g in \/gE = Since in the colour plate 
model, we set the zero of time not at maximum overlap, but at the instant when the plates 
have completely crossed each other, we shall take the initial energy density e = SOOGeV/ fim?. 
We shall also be guided by the results of Geiger and Kapusta in our choice of the values 
of Tc, and take = 0.2 fm, to be a realistic value The scale for the hydrodynamic 



limit will be set by the choice for Tc as well as the formation time tq. A typical value for the 



analysis at hand is 0.001 fm, which we employ here. Finally, we study the other extreme 
case, the coUisionless limit, by the choice Tc = 5/m. It is instructive to compare how the 
realistic regime behaves in relation to the results that we obtain for the two extreme limits. 

We have studied the decay of the chromoelectric field, the evolution of the particle 
energy and number densities, the evolution of the distribution function and its approach to 
the equlibrium state. We have compared them with the corresponding Abelian results. In 
addition, we have also evaluated a quintessentially non- Abelian quantity, viz, the expectation 
value of the angle that the quark charge makes with the direction of the chromoelectric field. 
In the current model the chromoelectric field not only decays to produce the qq pairs, but 
is itself built up by the receding nuclei which act as colour plates. We have, therefore, 
also studied the particle energy per unit transverse area as a function of ordinary time to 
highlight this feature. 



A. Discussion of the results 

We shall present the results of our analysis in the three regimes corresponding to the 
hydrodynamic, the realistic, and the coUisionless cases. Comparison between the Abelian 
and the non- Abelian systems will be taken up subsequently. 

Consider the hydrodynamic limit first. The physical quantity of utmost importance 
is the chromoelectric field, whose dynamics is most readily determined in the Flux tube 
model at hand. Indeed, apart from determining the production and the acceleration of the 
partons in the Flux tube model at hand, it has an additional role which has been emphasized 
by Svetitsky |^|: the charms which are produced in the pre-equilibrium stage would not 



only interact with the gas of light quarks and gluons, before either forming a J/ip or an 
open charm mesonic state, but will also be influenced by the background field. Thus the 
study the evolution of the mean electric field in the pre-equilibrium stage acquires an added 
importance. It is a merit of the Flux tube model that we can readily determine the evolution 
of the background electric field. Note that in contrast, no such information can be extracted 



in the more microscopic models such as that of Geiger and Kapusta |T^. In fig. 1 the decay 



of the field is shown. Recall that Tc = .001/m. As can be seen from Fig. 1, the electric field 
has hardly decayed at all, with a percentage decay less than 2 %, even at r = 1.5/m. The 
corresponding particle energy density, shown in Fig. 4 is also negligibly small, with the ratio 
with respect to the initial energy density being ~ 10^^. The Abelian situation is an order 
of magnitude better, but is still hopelessly small. And indeed, the number density is also 
very small, and as shown in Fig. 7 gets saturated at ~ .05//m^. and the temperature also 
stabilizes to a value ~ AOMeV, which is much less than the temperatures required. Clearly, 
it is very unlikely that the plasma does not go through a prequilibrium phase. 

The behaviour of the system for Tc = 0.2 and 5/m is in sharp contrast to the case 
discussed above. There is a significant decay of the electric field, which is shown in Figs. 
2 and 3 for the two respective values of Tc. The electric field decays by about 15% of the 
original value, at r = 1.5/m. The corresponding energy densities are also not very different, 
and yield ~ 10 — 15% of the initial field energy density (see Figs. 5 and 6). The corresponding 
temperatures are ~ SOOMeV, which is quite realistic. Note that the parton cascade model 
||T7| predicts a value about 300 MeV at 2. 5/m. The real difference between the collisionless 
case and the "realistic" case is in the number density. Indeed as Figs. 8 and 9 show, the 
plasma produced for Tc = 0.2 fm has a number density ~ 15 — 20//m^, the corresponding 
number for Tc = 5fm is three to four times smaller. In other words, the system is more dense 
in the former case than the latter. Clearly, this distinction should show up in signatures 
such as dilepton production which are sensitive to both the number as well as the energy 
density. 

It thus follows that the flux tube model, with the incorpration of the colour degrees 
of freedom, deflnitely rules out instantaneous equilibration as envisaged in the approach 
of Baym 0; further, it also distinguishes the collisionless limit from realistic values of 
equilibration time. 

Finally, it should be emphasized that the flux tube model as employed here not only 
pumps in the fleld energy to particle production, but also contributes to the fleld energy 



by virtue of the recession of the colour charged nuclei from each other. In fact, a crude 
estimate of the time required to convert all the plate energy to the field energy turns out 
to be ~ 5/m for 200G'eV^/nucleon. In any case, it is therefore misleading to interpret the 
ratios we have shown in Figs. 4-6 as the fraction of the total field energy that has gone into 
the particles. To emphasize this, we have evaluated the dependence of the field and particle 
energies per unit area by integrating over the contribution along the longitudinal direction. 
The results are shown in Fig. 10, from which it is clear that there is a lot more energy in 
the field than in the particles. For the same reason, it is also misleading to interpret the 
field energy density at r = as the counterpart of the energy densities employed in lattice 
computations to study the transition from the hadronc to the QGP phase. 



B. Comparison with Abelian results 

It is clear from Figs. 1-9 that the Abelian and the non-Abelian results bear little re- 
semblence, belying the expectation that the "Abelian dominance" holds in this case. The 
difference is most dramatically highlighted at = 0.2/m. Whereas almost all the initial field 
has decayed in the Abelian case, only 30% has done so in the non-Abelian case. Accordingly, 
The particle energy density is larger by a factor of ~ 5 , the number density by a factor of 
~ 4, and yields an abnormally large value of~ SOOMeV for the temperature. Earlier Abelian 
calculations |[T^JT5|JT6[] yielded reasonable values for the temperature because of unrealistic 



values for the initial field energy. It is a universal feature that the non-Abelian plasma is 
rarer and cooler than its Abelian counterpart. There are other significant features. Consider 
the hydrodynamic case. It is seen that although the decay in the non-Abelian field is more 
than the Abelian field, the corresponding energy density is smaller than in Abelian case. 
Indeed, in the hydrodynamic limit, the Abelian ananlysis yields a temperature ~ 200MeV, 
in contrast to ~ AOMeV in the non-Abelian case. Thus the incorporation of the colour 
degree of freedom rules out instantaneous hydrodynamic evolution. Of course, a colourless 
plasma does not have reasonable temperature for any other Tc- 



Yet another interesting aspect that emerges from our studies is in the close interplay 
between the value of Tc and the colour degree of freedom. While the field decays faster in 
the non-Abelian case in the hydrodynamic limit, the trend reverses at Tc = 0.2/m and gets 
restored in the coUisionless limit. In contrast, the particle energy density is smaller for a 
coloured plasma in hydrodynamic limit. Further, it continues to be so at Tc = 0.2/m but 
becomes larger than the Abelian case in the coUisionless limit. These non-trivial manifesta- 
tions of the colour charge may be expected to have an important bearing on the other bulk 
properties of the plasma. 

In order to gain some insight into the above features, we have also evaluated the expecta- 
tion value of © =< cos^O > which yields the rms value that the particle charge makes with 
the field direction in the colour space. Note that the above quantity is gauge invariant, and 
hence physical. Fig. 11 shows 6(t) for Tc — 0.2 fm. It may be seen that the value saturates 
around 0.25, corresponding to ^ ~ 7r/3. The corresponding value in the non-Abelian case 
is strictly zero. The effect of the background field is thus reduced, leading to dominance of 
the collision term in the expansion of the particles. Indeed, it shows up most clearly in the 
approach to the equilibrium, where one now expects that the plasma equilibrates the fastest 
in the direction normal to that of the field (in the colour space). This is corroborated as 
may be seen in Fig. 12. 

V. CONCLUSION 

To conclude, we have studied the production and the equilibration of a genuinely non- 
Abelian plasma with the colour degree of freedom incorporated in both the source and the 
background field term in the transport equation. In the SU (2) gauge theory that we have 
considered here, the distribution function is defined in the extended phase space. We find 
that this approach recaptures in an elegant manner many of the findings of the more micro- 
scopic parton cascade model. It has the further advantage that it indeed exhibits the colour 
degree of freedom manifestly, and allows us to compute various gauge invariant quantities. 
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Significantly, we also find that the Abelian approximation, employed hitherto in most stud- 
ies of equilibration of QGP is rather too drastic to be used for any quantitative analysis and 
comparison with the experiments. The study also almost rules out instantaneous equilibra- 
tion, and also strongly suggests that the coUisionless hmit may also not be the favoured in 
URHIC. 

To be sure, we have not made any comparison with the experimental findings here, for 
the simple reason that we are as yet dealing with a simpler gauge group , and ignored the 
gluonic component altogether. The indications from the present study arc unmistakable, 
though. Indeed, the particle production is enhanced because of increase in the phase space 
available, and for the same reason the plasma will be cooler than the Abelian counterpart. 
The energy goes to the colour degree of freedom, and docs not simply heat the system as 
it would happen in a colourless plasma. If we consider the realistic SU (3) case, this feature 
will get further accentuated; For the same initial configuration, we may expect a rarer and 
a cooler plasma. Of course, there are other features which are intrinsic to SU {?>)■. there is 
yet another Casimir invariant, and the gluon term will also have to be incorporated. These 
are under study at the present and will be reported else where, with a full discussion of the 
signatures of the QGP in the flux tube model. 
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Figure captions 



FIG. 1. Decay of the chromoelectric field as a function of proper time (in units of fermi), 
for Tc = .001/m. The sohd hne refers to the non-Abehan case, and the broken hne to the 
Abehan case. 

FIG. 2. Decay of the chromoelectric field as a function of proper time (in units of fermi), 
for Tc = .2fm. The solid line refers to the non-Abelian case, and the broken line to the 
Abehan case. 

FIG. 3. Decay of the chromoelectric field as a function of proper time (in units of fermi), 
for Tc — 5fm. The solid line refers to the non-Abelian case, and the broken line to the 
Abehan case. 

FIG. 4. The particle energy density scaled w.r.t the initial field energy density as a function 
of proper time (in units of fermi), for Tc = .001/m. The solid line refers to the non-Abelian 
case, and the broken line to the Abelian case. 

FIG. 5. The particle energy density scaled w.r.t the initial field energy density as a function 
of proper time (in units of fermi), for Tc = .2/m. The solid line refers to the non-Abelian 
case, and the broken line to the Abelian case. 

FIG. 6. The particle energy density scaled w.r.t the initial field energy density as a function 
of proper time (in units of fermi), for Tc = 5fm. The solid line refers to the non-Abelian 
case, and the broken line to the Abelian case. 

FIG. 7. The particle number density as a function of proper time (in units of fermi), for 
Tc = .001/m. The solid line refers to the non-Abelian case, and the broken line to the 
Abehan case. 

FIG. 8. The particle number density as a function of proper time (in units of fermi), for 
Tc = .2/m. The solid line refers to the non-Abelian case, and the broken line to the Abelian 
case. 

FIG. 9. The particle number density as a function of proper time (in units of fermi), for 
Tc = 5fm. The solid line refers to the non-Abelian case, and the broken line to the Abelian 



case. 

FIG. 10. The particle energy/unit transverse area (solid line) and the field energy/unit 
transverse area (dashed line) as a function of time (in fermi). The energy densities are in 
GeV/fw? and = .2/m. 

FIG. 11. < cos^O > as a function of proper time (in fermi) at = .2/m. 
FIG. 12. f/feq as a function of proper time at pt — 200MeV, ^ — and Tc = 0.2 fm for 
three different angles corresponding to cos6 — (solid line), cos9 — .25 (dash fine just below 
the solid line), and cos6 = 1 (the other dash line). Note that the equlibration is fastest at 
e = 7r/2. 



